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Abstract We study a thermodynamically consistent implementation 
of the nucleon volume in the mean field theory, and find that this volume has 
large consequences on the properties of hadronic matter under extreme con- 
ditions such as in astrophysical objects and high energy heavy-ion collisions. 
It is shown that we can reproduce the critical temperature T c ~ 200 MeV 
predicted by lattice QCD calculations for the phase transition from hadronic 
matter to quark-gluon plasma, by using parameters which are adjusted to 
fit all empirical data for normal nuclear matter. 
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One of the main goals of nuclear physics is to study the properties of 
hadronic matter under extreme conditions, such as in high density and high 
temperature astrophysical objects and high energy heavy-ion collisions [|J|, 
based on our knowledge of normal nuclear matter properties. There are some 
difficulties to reproduce simultaneously all known properties for nuclear mat- 
ter under normal and extreme conditions by using the naive conventional 
scenario of point-like hadrons. It has been observed § that in the framework 
of the local relativistic mean field theory the mass limit for neutron stars is 
too low when using a nuclear equation of state with a realistic compression 
modulus K around 200 MeV. Also, there is no reasonable phase transition 
to the quark-gluon plasma for vanishing net baryon number at large T due 
to the possible excitation of a vast number of hadronic resonances H, |j. The 
purpose of this letter is to show that the above difficulties can be removed 
simultaneously by simply taking into account the effect from the finite vol- 
umes of nucleons. That such a consideration is reasonable can be motivated 
by the success of the Van der Waals equation of state to describe proper- 
ties of real gases. The limits of the validity of the effective field theory of 
point-like mesons and baryons are also indicated by the off-shell behaviors of 
hadronic structure functions due to the internal quark and gluon structure of 
hadrons ||. By adopting a thermodynamically consistent procedure in the 
mean field theory, we have successfully produced § a mass limit for neutron 
stars which is compatible with the observation, while still kept the compres- 
sion modulus for normal nuclear matter around its empirical value K = 240 
MeV. In this letter we will focus our attention on the phase transition to 
the quark-gluon plasma which is expected to occur in high energy heavy-ion 
collisions 0. We will show that we can reproduce the critical temperature 
T c ~ 200 MeV predicted by lattice calculations of pure SU(3) gauge theory 
[]l], ||, H using the parameters which are adjusted to fit all empirical data for 
normal nuclear matter |i~0|] , namely, the equilibrium density po = (4^r 3 ) -1 
with r = 1.175 fm, the binding energy per nucleon B/A = 15.986 MeV, the 
asymmetry energy E as = 36.5 MeV, and the compression modulus K = 240 
MeV, with the bag constant for the quark-gluon plasma phase ranging from 



2 



B 1 / 4 = 125 MeV to about 300 MeV. 

In our approach hadronic matter and the quark-gluon plasma are de- 
scribed separately by different equations of state and a transition is obtained 
via Gibbs' phase coexistence conditions. The effect of the finite volume has 
been extensively discussed in the ideal hadronic gas model [11] which cannot 
reproduce the nuclear ground state properties. We introduce the finite vol- 
ume of nucleons in a thermodynamically self-consistent way j|, ||, [l^, 13] in 
the mean field theory framework, which has been widely used for the descrip- 
tion of nuclear matter, finite nuclei, and nuclear dynamics [[D|]. One new 
aspect of our approach is that the equation of state for the hadronic phase 
is applicable to hadronic matter under extreme conditions as well as normal 
conditions, although there are problems with causality at very high density 
[|], [|. The quark-gluon plasma is treated as a gas of relativistic massless 
quarks, antiquarks, and gluons, including first-order QCD corrections and 
non-perturbative effects by an overall bag constant B. 

For the sake of simplicity, we assume symmetric, homogeneous, and 
isotropic hadronic matter with Nb baryon species (with equal properties) in 
a volume V, and restrict ourselves to non-strange degrees of freedom. In the 
ideal gas model the pressures of nucleons and antinucleons can be written 
as 

kT 

V° N (^T) = — CGu.T); (1) 
kT 

VUn,T) = — C(-M,T), (2) 
with T) defined by 



F -C(M,T) =N B J j^ln[l + e^-% (3) 

where j3 = 1/kT with T being the temperature, /x is the chemical potential 
and e = \Jp 2 + x 2 is the energy of a nucleon of 3-momentum p and mass \. 
We now consider a system of nucleons and antinucleons with finite volume 
r. Assuming that the repulsive force due to the finite volume is operative 
between nucleons or antinucleons but not between a nucleon and an antinu- 
cleon due to the possibility of matter-antimatter annihilation, we can simply 
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discuss the finite volume effect for nucleons and antinucleons separately. The 
pressures can be written as Q [l3|] 

kT 

V n (jm,T) = V%{P,T) = — C(A,T); (4) 



V^,T)=V%(H',T) = — a-fi',T), (5) 



kT 
~V 

with the effective chemical potentials p for nucleons and p' for antinucleons 
given by 

p, = fj, - tVn(h,T); p' = n + TVfj{p,,T), (6) 

as required to keep thermodynamical consistency |l^] . In the mean field 
theory of Bodmer and Boguta |l4|], the energy density of the hadronic phase 
is 



m 
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£h = £n + £fi + U + 2ita UJ ^TP h ; ( 7 ) 

with U given by 



U = -^(1 - x) 2 [l + a 1 (l- X ) + « 2 (1 - X?l (8) 
direct 

where ph is the net baryon number density, m and m w are the masses of 
free nucleon and u meson, and a, ai, 02, and ot w are the mean field theory 
parameters adjusted by fitting the normal ground state properties with the 
finite volume of the nucleon r = ^-o? (with a = 0.62 fm) ||. The effec- 
tive chemical potentials v for nucleons and v' for antinucleons are further 
modified to 

2 2 

Try 7 n 

v = p - A^oiu^-ph] v = p' - 'iira^^-ph (9) 



due to the effects from mean fields |16|]. The baryon number and kinetic 
energy densities for nucleons and antinucleons are 

_ r 2N B d 3 p 1 - PN r 

PN J (2vr)3 e P{y/M-v) + i [ ' 

= f 2N B d 3 p 1 - Pn t 
Pn J (2vr)3 e /3(V^?+,') + i' 1 ' 
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„ _ , 2N B d 3 pV^Tx I (l-pNT) 
N ~ 1 {2nf ePiy/M-v) + i ' 1 ' 



N J (2vr)3 e/3 (V^?+,') + 1 ' [6> 
The kinetic pressures for nucleons and antinucleons can be written as 

hT 

V N ^T)=r%{u,T) = —Ciy,T)- (14) 

kT 

V*(p,T)=V%(v',T) = — C(-u',T), (15) 

Thus we obtain the total pressure and the net baryon number density for 
the hadronic phase: 



in 
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V h {n,T) =p N (ji,T)+Vjf(ii,T)-U + 2ica u — I ti; (16) 

Ph = PN-PR- ( 17 ) 

The effective nucleon mass x is fixed by minimizing the energy density at 
constant baryon number density p^ and entropy density; which is equivalent 
to maximize the pressure at fixed T and p. For the quark-gluon plasma 
phase, the pressure, energy density and net baryon number density are given 
by § 

-r> r rp \ Svr 2 ^ 4 15a s 7n%* 50a, /#T| p\ 

V qgp {p q ,T q ) = _(1 - _) + ^/[-go-(l - + (— + ^) 

2rv 

)}~B; (18) 

7T 

Svr 2 ^ 15a 77r 2 T 4 5Qa ujlj „} 2 

" -T5- (1 -^ )+iV/[ ^0- (1 -n7 )+3( ^ + 4^ )(1 - — )]+ ^' 

(19) 

^ = ^W 2 + 4)(1-^), (20) 
where B is the bag constant and a s is the QCD running coupling constant 
which depends on the quark-gluon plasma temperature T q and the quark 
chemical potential p q through [17] 

4-7T 

a s (p q ,Tg) = (n _ 2 N f /3)ln[(0.8p 2 + 15.622T 2 )/A 2 ] ' (21) 
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with A being the QCD scale parameter. 

To examine the phase transition between hadronic matter and the quark- 
gluon plasma, we apply Gibbs' conditions of phase coexistence, i.e., Vh = 
Vqgp] T = T q ; and p = 3p q . As a first approximation, we calculate the 
simplest case with Nb = 2 (i.e., only protons and neutrons) and Nf = 
2 (i.e., only u- and d- quarks). The range of B is chosen from B 1 ^ = 
125 MeV to about 300 MeV which is consistent with the results from a 



bag model analysis of hadron spectroscopy [g, 11, 15, IS]. The QCD scale 
parameter is now quite well determined to be A ~ 200 MeV |20]]. Fig. 1 
presents the numerical result for the T — p phase diagram with B 1 ^ = 200 
MeV. In Fig. 2, the dependence of the net baryon number densities on 
the temperature for the hadronic matter and the quark-gluon plasma phase 
boundaries are plotted. We see that we can reproduce a reasonable phase 
transition temperature of T c = 170 MeV for vanishing net baryon number 
density. T c changes from about 120 to 240 MeV by changing the parameter 
B V4 f rom 12 5 MeV to 300 MeV. Changing a = 0.62 fm to a = 0.74 fm (and 
readjusting the mean field parameters) causes a very small change (within 
1 MeV) in T c . By turning off the 0(a s ) correction we find T c = 140 MeV, 
thus the first order perturbative QCD corrections change the ideal quark- 
gluon gas result for the phase transition temperature by 20%. The finite 
volume effect is very large for nucleons and negligible for antinucleons, as 
can be seen from Fig. 1 by comparing the effective chemical potentials v for 
nucleons and v' for antinucleons with the effective chemical potential 

t 

,2 



TfT 

(jl- Anauj—^ph, (22) 



where the contribution from the finite volume is absent. We see that v 
is reduced by about 40% in comparison with v in the high density region 
whereas v 1 coincides everywhere with v. 

However, the finite volume effect should also play a non-trivial role for 
antibaryons if there is a large number of hadronic resonance excitations. It 
has been indicated ||, H that in a point-like hadronic scenario the pressure of 
hadronic matter will be larger than that of the quark-gluon plasma at large 
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T and vanishing ph due to the large number of hadronic resonances, thus one 
cannot reproduce a phase transition to the quark-gluon plasma as predicted 
by lattice QCD calculations. In our approach this situation will not occur 
due to the finite volume of baryons and antibaryons. To reflect this aspect 
of the approach, in a first rough estimate we simply use a large Nb, which is 
equivalent to assume a large number of baryons with the same properties as 
nucleons. We find that a phase transition still occurs with T c ~ 260 MeV, 
even if we use a very large Nb of 2 7 = 128. This means that by taking into 
account the finite volume of baryons we still have a transition to the quark- 
gluon plasma even if there are a large number of hadronic resonances, thus we 
overcome the difficulty of the point-like treatment of hadronic matter. From 
the T — [i phase diagram (Fig. 3), where the effective chemical potentials 
for baryons v and antibaryons v' are also plotted, we see that the finite 
volume effect is also large for antibaryons at large T, in contrast with the 
simple case of only considering nucleons and antinucleons. The necessity 
of considering other baryonic degrees of freedom can be also seen from the 
large difference between the two T — fj, diagrams Figs. 1 and 3: the effective 
chemical potentials for baryons v is much more reduced in the latter case 
than in the former. 

From Fig. 2 we see that there are large differences between the baryon 
number densities for the hadronic matter and the quark-gluon plasma at 
the phase transition at given T. This aspect does not change in the case 
of a large baryon number Nb = 128, as can be seen from Fig. 4. Another 
interesting feature of our results is that the effective mass of the nucleon 
is non-zero in the whole hadronic region: we find x — 0.6m for the case of 
Fig. 1 and x — 0.5m for Fig. 3. This is different from previous expectations 
of a strong reduction of the effective baryon mass near T c [ 21 1 but consistent 
with the results of ref. Q. Of course, the exact quantitative details depend 
on the specific couplings of other baryonic resonances to the a field and 
their effective masses. Such detailed investigations are out of the scope of 
the present letter. 

In summary, we have studied the finite volume effect of nucleons by 
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using a thermodynamically consistent procedure in the mean field theory 
approach, and found that this volume effect has large consequences in high 
density and high temperature nuclear matter concerning the limit of neu- 
tron star masses and the quark-gluon plasma phase transition in high energy 
heavy-ion collisions. Further studies should include a realistic mass spec- 
trum of baryonic resonances, instead of our simple estimate with 2 7 equal 
resonance species, and extend the discussion to strange degrees of freedom. 

We would like to acknowledge many helpful discussions with Jiirgen 
Schaffner, Debades Bandyopadhyay, Raffael Mattiello, and Horst Stocker. 
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Figure Captions 

Fig. 1. The T — \i diagram for the phase equilibrium with A = 200 MeV and 
B 1 / 4 = 200 MeV. The solid curve is the chemical potential fi. The 
dashed and dot-dashed curves are the effective chemical potentials v 
for nucleons and v obtained by turning off finite volume corrections 
(see eq.(24)). The effective chemical potential v' for antinucleons is in 
coincidence with v. 

Fig. 2. The solid and dashed curves are the baryon number densities ph and 
p q gp for hadronic matter and the quark-gluon plasma under phase 
equilibrium with the same parameters as in Fig. 1. 

Fig. 3. Same as Fig. 1, but with a very large Nb = 128. The dotted curve is 
the effective chemical potential v' for antibaryons. 

Fig. 4. Same as Fig. 2, but with the same parameters as in Fig. 3. 
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